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SOFTWARE  FOR  A  COMPUTER  CONTROLLED  LASER  DOPPLER  VELOCIMETER 


I.  INTRODUCTION 

The  hazard  associated  with  fire  in  confined  spaces  is 
critical  to  the  Navy  due  to  the  vulnerability  of  a  ship  to  the 
disastrous  effects  of  an  uncontrollable  fire.  An  understand¬ 
ing  of  the  physical  mixing  processes  in  diffusion  flames, 
which  are  characteristic  of  most  unwanted  fires,  is  neces¬ 
sary.  By  knowing  the  mixing  patterns,  a  better  understanding 
of  fire  growth  in  confined  areas  can  be  developed.  This 
ultimately  will  lead  to  ways  of  controlling  materials, 
geometry,  design  and  structure  to  fire  harden  our  ships  more. 
The  need  for  better  mixing  information  has  prompted  research 
efforts  toward  the  development  of  a  technique  to  measure 
turbulence  in  combustion  processes  using  laser  doppler 
velocimetry  (LDV) . 

The  LDV  is  a  state-of-the-art  diagnostic  tool  to  measure 
the  velocity  of  particles  moving  in  a  flow  field.  The  ability 
to  characterize  highly  turbulent  flow  fluctuations  at  near 
zero  mean  velocity  and  to  make  measurements  without  obstruct¬ 
ing  the  flow  makes  turbulence  measurement  via  LDV  preferable 
to  that  of  other  developed  methods,  such  as  the  hot  wire 
anemometer.  The  inability  of  the  hot  wire  anemometer  to 
measure  turbulence  in  smoke-laden  environments  makes  the  LDV 
approach  particularly  attractive  (1,2)  . 

The  principle  behind  laser  doppler  velocimetry  lies  in 
the  phenomenon  that  particles  suspended  in  a  flow  field  will 
scatter  laser  light  that  is  directed  upon  them.  The  frequency 
of  the  scattered  light  depends  on  the  particle  velocity  and 
the  angle  through  which  the  light  has  been  scattered.  The 
frequency  of  the  scattered  light  is  different  from  the  fre¬ 
quency  of  the  incident  laser  light.  This  frequency  difference 
is  the  doppler  frequency,  which  can  be  detected  by  electronic 
processing. 

This  is  a  continuing  research  effort,  the  result  of 
which,  to  date,  has  been  the  development  of  computer  software 
to  control  the  LDV  system.  An  example  run  follows,  for 
illustrative  purposes  in  demonstration  of  the  software 
capabilities. 

Manuacript  approved  March  22,  1984. 
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II.  EQUIPMENT 


The  LDV  System  (Disa  Electronics)  is  comprised  of  a  two- 
color  laser  and  associated  optics ,  two  counter  processors , 
electronic  frequency  shifters/  and  a  buffer  interface.  Figure 
1  illustrates  the  overall  schematic  of  the  LDV  system.  The 
pieces  of  the  equipment  describing  the  LDV  System  are  des¬ 
cribed  below. 

II.  A.  LASER 

The  blue/green  colored  beam  from  a  four  wattf  argon  ion 
laser  is  passed  through  a  series  of  color-neutral  and  color- 
separating  beam  splitters  to  create  three  beams,  blue/green, 
green  and  blue,  as  shown  in  Figure  2.  The  three  beams  are 
focused  by  the  transmitting  optics  to  an  intersection  point, 
approximately  16  feet  from  the  laser  where  the  fluid  whose 
flow  patterns  are  being  measured  is  located.  The  purpose  of 
choosing  such  an  unusually  large  focal  length  (most  LDV 
systems  employ  a  three  to  four  foot  focal  length)  is  to  have 
the  capability,  when  necessary,  to  make  measurements  during 
large  scale  fire  tests  where  it  is  desirable  to  keep  the 
hardware  at  a  sufficiently  large  distance  from  the  fire 
source. 

The  intersection  of  the  green  and  blue/green  beams 
results  in  the  formation  of  interference  fringes.  A  particle 
from  the  flow  field  which  passes  through  the  interference 
fringes  scatters  light,  the  intensity  of  which  rises  and  falls 
as  the  particle  moves  from  bright  fringe  to  dark  fringe.  This 
scattered  light  is  collected  on  a  photomultiplier  tube  in  the 
backscatter  mode  and  converted  to  a  frequency-modulated  cur¬ 
rent  burst.  Further  processing  occurs  in  the  counter  proces¬ 
sors  (as  discussed  below)  to  give  the  vertical  component  of 
velocity  of  the  particle  in  the  flow  field.  The  intersection 
of  the  blue  and  blue/green  beams  similarly  creates  interfer¬ 
ence  fringes  from  which  a  particle  passing  through  scatters 
light.  This  scattered  light  is  collected  on  a  second  photo¬ 
multiplier  tube,  converted  to  a  doppler  signal,  and  further 
processed  in  the  second  counter  processor  to  give  the 
horizontal  velocity  component  of  the  same  particle  in  the  flow 
field. 

II.  B.  COUNTER  PROCESSOR 

The  frequency-modulated  current  burst  from  each  of  the 
photomultipliers  is  amplified  and  band-pass  filtered  in 
separate  counter  processors.  Each  burst  represents  the 
passage  of  a  single  particle  and  contains  information  on  the 
number  of  fringes  passed  by  the  particle.  The  frequency  of 
the  doppler  burst  signal  is  detected  by  internal  logic 
circuitry.  The  distance  between  the  fringes  is  known  from  the 
geometry  of  the  optical  beams  and  the  wavelength  of  the  laser 
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line.  Instantaneous  particle  velocity  may  then  be  calculated 
by  multiplying  the  fringe  spacing  by  the  measured  doppler 
signal  frequency. 

The  counter  processors  have  the  capability  to  send  four 
types  of  data  to  the  buffer  interface.  Doppler  frequency  data 
and  sample  interval  time  (time  interval  between  each  succes¬ 
sive  validated  signal)  data  are  contained  on  the  pins  of  one 
of  the  rear  panel  multi-connectors  of  the  counter  processor. 
The  other  gives  information  on  the  number  of  fringes  and  burst 
time,  i.e.,  time  between  each  doppler  burst. 

II.  C.  BUFFER  INTERFACE 

The  Buffer  Interface  acts  as  a  data  and  control  interface 
in  the  transfer  of  digital  data  between  the  counter  processors 
and  the  Hewlett-Packard-1000  mini-computer.  It  functions  as  a 
buffer  memory  for  the  randomly  arriving  data  from  the  counter 
processors  before  it  is  transferred  to  the  HP-1000.  It  also 
receives  the  control  words  from  the  computer,  as  written  in 
the  software,  so  that  the  LDV  will  operate  with  the  appropri¬ 
ate  data  and  control  channels  opened. 

Synchronization  within  a  selected  time  interval  between 
the  data  from  the  two  counter  processors  is  achieved  with  the 
coincidence  filter  board  inside  the  Buffer  Interface.  The 
selection  of  input  and  output  channels  and  the  setting  of 
coincidence  time  intervals  are  all  software-controlled. 

The  Buffer  Interface  has  been  designed  to  interface  with 
the  DEC  PDP-11  mini-computer  family.  An  additional  circuit 
board  was  added  to  the  Buffer  Interface  which  converts  the 
control  and  data  format  to  that  of  the  HP  Microcircuit  Board 
12566B  for  Hewlett  Packard  mini-computer  in  the  2100  series. 

Data  are  transferred  from  the  additional  interface  board, 
with  a  16  bit,  parallel  full  duplex  format,  to  the  12566 
Board  of  the  HP  21MX  series  mini-computer.  The  12566  Board 
has  16  bit  input  and  output  registers.  A  system  of  electronic 
"handshakes"  controls  the  signal  and  data  lines,  as  described 
in  the  software  section  of  this  report. 

III.  SOFTWARE 

The  transfer  of  information  to  and  from  the  LDV  is  soft¬ 
ware  controlled.  Control  of  the  LDV  by  the  HP-1000  computer 
is  achieved  by  a  software  driver.  The  driver  is  called  by  the 
main  FORTRAN  program,  program  LDV,  both  of  which  were  written 
at  the  Naval  Research  Laboratory  (NRL)  and  are  discussed  in 
the  following  sections. 
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III.  A.  DRIVER 


A  relocatable  HP-assembly  language  driver  was  written  to 
enable  the  computer  to  communicate  with  the  LDV  Buffer 
Interface.  The  driver,  as  shown  in  Appendix  A,  is  FORTRAN 
callable,  and  contains  all  of  the  necessary  "set  control," 

"set  flag,"  etc.  instructions  to  control  and  read  data  from 
the  LDV. 

A  non-interrupt  data  transfer  method  was  chosen  due  to 
the  higher  degree  of  ease  with  which  this  type  of  driver  could 
be  written  compared  to  that  of  an  interrupt-type  driver.  The 
non-interrupt  method  involves  a  "wait-for-flag"  routine  when 
inputting  LDV  data  into  the  data-register .  The  computer 
signals  the  interface  card  to  start  the  LDV  device.  The 
driver  loops  until  the  LDV  is  ready  to  input  data  into  the 
data-register.  Readiness  to  input  data  is  expressed  by 
setting  a  flag  flip-flop.  After  the  flag  is  set,  the  program 
continues.  At  times,  the  doppler  signal  is  of  insufficient 
strength  to  cause  the  flag  to  set.  To  prevent  the  driver  from 
looping  indefinitely,  a  routine  to  wait  1000  times  for  the 
flag  was  incorporated.  If  the  flag  does  not  set  after  1000 
attempts,  the  program  is  terminated.  The  wait-for-flag 
routine  is  not  required  when  outputting  control  words  from  the 
computer  to  the  LDV. 

The  memory  protect  logic  feature,  which  is  incorporated 
in  the  HP-1000  Real  Time  Executive  operating  system,  limits 
control  of  all  Input/Output  operations  to  interrupt  control 
only.  Thus,  the  memory  protect  fence  is  disabled  by  assembly 
language  software,  so  that  the  non-interrupt  driver  has 
control  of  the  LDV. 

The  driver  consists  of  two  main  segments,  the  initiation 
section  and  the  continuation  section.  The  initiation  section 
performs  two  basic  functions:  (a)  to  clear  the  buffer  memory 
before  a  new  record  of  data  is  to  be  entered,  and  (b)  to  send 
the  control  words  to  the  LDV. 

The  collection  of  data  is  accomplished  in  the  continua¬ 
tion  section.  Counting  variables  are  incremented  each  time  a 
data  word  is  collected  so  that  data  rate  can  be  acknowledged. 
The  memory  protect  fence  is  re-enabled  after  a  buffer  of  data 
is  collected.  This  gives  the  interrupt  system  control  of  the 
input/output  peripheral  devices  before  returning  to  main 
FORTRAN  program.  Three  control  words  are  sent  to  the  LDV 
which  open  up  appropriate  channels  for  data  collection.  The 
control  words,  which  are  collected  in  the  initiation  section 
of  the  driver  and  are  sent  out  in  the  main  program,  perform 
three  functions. 


The  first  word  turns  on  the  coincidence  filter  board  and 
opens  up  channels  1  and  2.  The  second  control  word  sets  the 
time  interval  between  successive  samples  to  125  microseconds. 
This  will  ensure  that  the  data  from  the  two  counter  processors 
will  be  synchronized  within  this  time  interval.  The  third 
control  word  is  the  static  control  word  and  must  remain  set 
during  operation.  Channels  1  and  2  (connected  with  the 
counter  processors)  and  channel  4  (connected  with  the 
Coincidence  Filter  Board)  are  opened  for  input.  Doppler 
frequency  data  and  sample  interval  data  may  be  simultaneously 
collected  from  one  of  the  multi-connectors  on  the  counter 
processor;  fringe  number  and  burst  time  data  are 
simultaneously  collected  from  the  other. 

Because  the  12566  Board  does  not  give  a  data  transmitted 
signal  after  a  data  word  has  been  read  into  the  computer,  this 
signal  is  created  artificially  by  transmitting  a  device 
command  signal  with  the  address  of  octal  60000.  In  both  the 
initiation  and  continuation  sections,  the  number  60000  (octal) 
is  loaded  into  the  B  register.  The  instructions  OTB  LDV 
(output  the  contents  of  the  B  register  to  the  LDV)  and  STC 
LDV  (set  the  control  flip  flop  on  the  LDV  interface  card  in 
the  computer  I/O  cage)  result  in  the  artificial  creation  of  a 
data  transmitted  signal. 

III.  B.  PROGRAM  LDV 

The  main  FORTRAN  program,  as  shown  in  Appendix  B,  is 
comprised  of  three  main  sections  which  contain  calls  to  the 
assembly  language  driver  and  calls  to  several  subroutines. 

The  first  segment  is  used  to  test  the  LDV  by  taking  samples 
and  writing  this  out  to  any  terminal.  The  user  may  specify 
any  control  words  that  are  desired.  This  is  a  good  way  to 
test  to  see  if  the  different  channels  which  are  opened  are 
actually  sending  data. 

The  second  segment  is  used  to  collect  data  for  a  user- 
specified  amount  of  time  and  writes  the  raw  data  to  a  created 
disc  file.  The  maximum  number  of  512  data  transfer  words  can 
be  collected  in  a  record,  as  this  is  the  size  of  buffer  memory 
in  the  Buffer  Interface.  The  number  of  data  words  collected, 
represented  by  the  variable  ICNT,  is  one  of  the  variables 
passed  to  the  driver  from  the  main  program,  which  enables  the 
calculation  of  data  rate. 

Segment  number  three  is  used  to  sort  and  assemble  the  raw 
data  collected  on  the  disc  file  created  by  segment  number 
two.  A  second  data  file  is  created  which  contains  the  sorted 
data.  This  segment  reads  in  a  record  of  data  of  up  to  512 
words  from  the  raw  data  file.  Subroutine  ASMBL,  which  is  then 
called,  evaluates  each  data  word  based  on  the  type  of 
information  that  corresponds  to  each  of  the  16  bits  comprising 
the  data  words. 
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This  information  is  passed  to  subroutine  SORT  which 
identifies  the  data  as  either  doppler  frequency  data,  sample 
interval  data,  number  of  fringes  data,  or  burst  time  data. 
Doppler  frequency  data  is  converted  to  velocity,  by  the 
following  relationship: 

Velocity  =  Doppler  Frequency  x  Fringe  Spacing 


where 

x 

Fringe  Spacing  =  _ 

2  sin  9/2 

x  *  laser  wavelength 

e  =  angle  of  intersection  of  laser  beams 


Summations  of  the  velocity  data  are  performed  in 
Subroutine  STAT  for  the  purpose  of  calculating  average 
velocity  at  a  subsequent  point  in  the  program.  Subroutine 
GATHR  groups  the  velocity  and  sample  interval  time  data  by 
burst  and  writes  the  assembled  data  to  the  output  file.  At 
this  point,  program  statement  number  310,  the  program  loops 
back  up  to  program  statement  number  309  to  read  the  next 
record  of  data  of  up  to  512  words  from  the  raw  data  file. 

After  the  final  record  of  data  has  been  written  to  the 
output  file.  Subroutine  STAT  calculates  average  velocity  for 
each  of  the  counters,  i.e.,  horizontal  and  vertical  component, 
and  the  magnitude  and  direction  of  the  resultant  velocity 
vector. 

IV.  EXAMPLE  RUN 

The  following  example  is  intended  as  an  illustration  of 
the  utility  of  the  software  in  transferring  information  to  and 
from  the  LDV. 

The  vertical  and  horizontal  components  of  velocity  were 
measured  in  a  test  run  for  one  point  in  a  6  cm  methanol  pan 
fire.  The  laser  beam  intersected  at  a  vertical  position  1  cm 
above  the  lip  of  the  pan  and  in  the  center  of  the  pan  in  the 
horizontal  direction.  The  liquid  level  of  methanol  was  main¬ 
tained  constant  (3) .  Data  were  collected  for  approximately  30 
minutes. 

The  lack  of  formation  of  soot  particles  in  the  methanol 
flame,  due  to  methanol's  clean  burning  characteristics, 
creates  some  difficulty  with  LDV  measurement.  The  principle 
behind  the  LDV  technique  rests  in  the  fact  that  particles  in 
the  flow  field  will  scatter  laser  light  that  is  directed  upon 
them.  A  very  sparse  particle  concentration,  such  as  with 
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methanol,  results  in  a  very  low  data  rate  and  small  amplitude 
doppler  bursts.  A  doppler  signal  of  amplitude  less  than  200 
millivolts  will  not  be  validated  in  the  counter  processors, 
resulting  in  a  near  zero  data  rate.  The  unusually  long  focal 
length  of  16  feet  also  contributes  to  a  smaller  doppler  signal 
in  comparison  to  a  short  focal  length,  as  the  light  scattered 
from  the  particles  in  the  flow  field  travels  a  longer  distance 
before  being  collected  in  the  photomultiplier  tube.  The 
longer  distance  travelled  by  the  scattered  light  results  in 
reduced  signal  strength. 

In  this  run,  the  doppler  signal  was  boosted  by  adding 
small  amounts  of  n-heptane  to  the  fuel  pan.  This  resulted  in 
a  more  soot-laden  flame  of  orange  color.  The  flame  was  placed 
inside  a  wooden  box  with  several  small  cut-out  holes  for 
ventilation.  The  box  served  to  stabilize  the  flame  from  the 
air  currents  created  by  the  nearby  hood  and  air  ducts,  as  well 
as  to  elevate  the  smoke  particle  concentration  near  the  flame 
to  boost  the  doppler  signal  strength. 

Appendix  C  contains  a  listing  of  the  data  which  has  been 
grouped  by  doppler  burst.  Column  #2  lists  velocity  in  the 
vertical  direction  (meters  per  sec) ,  column  #3  lists  velocity 
in  the  horizontal  direction,  and  column  #4  gives  the  sample 
interval  time  for  that  doppler  burst  from  which  the  data  is 
generated.  For  the  vertical  component,  a  positive  number 
indicates  the  upward  direction,  while  negative  indicates 
downward.  A  positive  number  for  the  horizontal  component 
indicates  a  direction  to  the  right.  After  processing  all  of 
the  data,  the  unweighted  mean  velocity,  as  calculated  from  all 
the  data  points,  is  shown  along  with  the  angle  of  the  velocity 
vector,  in  reference  to  the  vertical  position. 

V.  CONCLUSIONS 

In  a  continuing  effort  to  understand  the  physical  mixing 
processes  in  diffusion  flames,  NRL  is  currently  involved  in 
the  development  of  a  technique  to  measure  the  velocity  profile 
in  combustion  processes  using  the  laser  doppler  velocimeter. 

To  date,  the  software  has  been  developed  for  a  computer- 
controlled  LDV.  Further  work  is  required  in  the  areas  of 
seeding  the  flow  field  to  boost  the  signal  to  noise  ratio  of 
the  doppler  signal  to  make  this  a  viable  research  tool. 
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APPENDIX  A  -  SOnVftKE  DRIVER 


J.DR’-'R 

T *0  y  uo  3 

IS  ON  CR00041  USING  00015  BLKS  R-0000 

0  40  1 

ASNB, R. L 

000  2 

HAH 

DR  VR  ,  7 

USER  HRITTEH  DR  I VER< SELECT  C0DE*10B) 

000  3 

ENT 

DR  VR 

0  00  4 

e:<t 

EMTR . 

8L1BR8LIBX 

o  oO  5 

* 

0006 

*  PINO  AOORESSES 

OF  PASSED  VARIABLES 

0  0«  7 

ICODE  NOP 

0  o  0  8 

IDaTA  NOP 

0  00  8 

I  £»  L  NOP 

0  010 

ISCU  NOP 

3011 

ISL  NOP 

0012 

ICHT  NOP 

0013 

IFIAC  NOP 

0014 

DPVR  NOP 

0015 

J  SB 

.  EHTR 

0016 

DEF 

ICQDE 

001? 

« 

0018 

*  DEFINE 

SELECT 

CODE  OF  12566  BOARD 

0  013 

LDV  EOU 

10B 

0  020 

* 

0021 

*  DETERMINE  TASK 

: 

0022 

LDA 

ICODE, 

1  I  CODE* 1  DO  INITIATION  SECTIOH  ONLY 

0423 

CPA 

TUO 

I  CODE* 2  DO  CONTINUATION  SECTION  ONLY 

0424 

JMP 

CONT 

I  CODE *3  00  BOTH  SECTIONS 

0025 

*******  ******************************  ****** I************************** 

0426 

•  INITIAT 

ION  SECTION: 

* 

0427 

********* 

******* 

0028 

LOB 

OC 

KEEP  60000B  IN  B 

0029 

LDA 

Kt 

SAVE  MEMORY  SIZE  TO  CLEAR 

0030 

CNA 

.  INA 

0031 

STA 

CNTl 

0032 

LDA 

ISL,  I 

USE  ISCU  LENGTH  AS  A  COUNTER 

0033 

CNA 

.  INA 

0034 

STA 

CNT2 

0033 

LDA 

ISCU 

SAVE  ADDRESS  OF  CONTROL  BUFFER 

0036 

STA 

PTR1 

0037 

* 

0038 

*  DISABLE 

THE  MEMORY  PROTECT  FENCE 

0039 

J  SB 

<L  IBR 

GO  PRIVILEGED 

0040 

NOP 

« 

0041 

* 

0442 

*  CLEAR  THE  USED 

BUFFER  MEMORY 

0043 

CLEAR  NOP 

0  044 

OTB 

LOV 

0045 

STC 

LDV.  C 

0046 

CLC 

LDV 

0047 

ISZ 

CNTl 

0048 

JMP 

CLEAR 

0049 

* 

0054 

*  SEND  OUT  CONTROL  MOROS  TO  INTERFACE 

0051 

SEND  LDA 

PTRl  .  I 

0452 

OTA 

LDV 

o  053 

STC 

LDV,  C 

0054 

CLC 

LDV 

0055 

IS2 

PTRl 

INCREMENT  BUFFER  ADDRESSES 

0056 

ISZ 

CNT2 

AND  COUNTER 

0  05? 

JMP 

SEND 

0053 

* 

11 


0439 

• 

RETURN 

IF  ICO DE • 1 

00*0 

LDA 

ICODE / I 

00*1 

CPA 

ONE 

00*2 

<1 NP 

BYE 

00*3 

*4 

1  **  4 

00*4 

# 

CONTINUATION  SECTION  i 

* 

00*3 

*4 

4 

■  +**•< 

00*6 

COHT 

NOP 

00*7 

LOA 

ONE 

ASSUME  A  SUCCESSFUL  CALL 

04** 

STA 

IFLAG • I 

00*9 

LDA 

IDL.  I 

USE  BUFFER  SIZE  AS  A  COUNTER 

0470 

CHA 

.  IHA 

0071 

STA 

CHT3 

0072 

LDA 

IDATA 

STORE  ADDRESS  OF  DATA  BUFFER 

0073 

STA 

PTR2 

0074 

LD8 

DC 

KEEP  60000B  IN  B-REG 

0073 

LDA 

ZERO 

INITIATE  DATA  COUNTER 

007* 

STA 

ICNT  >  I 

0077 

• 

0073 

« 

IF 

ICODE-3  HO  NEED 

TO  TURN  OFF  OPERATING  SYSTEM 

0079 

LDA 

ICODE / I 

0080 

CPA 

THRE 

0081 

4  NP 

MORE 

0082 

* 

0083 

4 

TURN  OFF  OPERATING 

SYSTEM 

0084 

4  SB 

SLIBR 

0083 

HOP 

0086 

* 

* 

0087 

* 

COLLECT 

OATA  TRANSFER  HORDS 

0488 

MORE 

HOP 

0489 

LDA 

THOU 

0490 

STA 

UA  IT 

0091 

OTB 

LDV 

0092 

STC 

LOV.C 

0093 

LIST 

HOP 

0094 

SFS 

LDV 

0093 

J  HP 

IMP 

0496 

LIA 

LDV 

0097 

CPA 

ZERO 

PREVENT  ZEROES  FROM  ENTERING 

0098 

•1  NP 

MORE 

BUFFER  MEMORY 

0099 

STA 

PTR2 »  I 

4100 

I  SZ 

PTR2 

0101 

HOP 

0142 

I  SZ 

ICHT,  I 

0103 

HOP 

0104 

ISZ 

CNT3 

0103 

4  NP 

MORE 

010* 

* 

0107 

• 

ENABLE 

THE  MEMORY 

PROTECT  FENCE  l  RETURN 

0108 

BYE 

HOP 

0109 

CLC 

LDV 

0110 

CLF 

LDV 

0111 

4  SB 

ILIBX 

GO  UNPRIVILEGED 

0112 

OEF 

*♦1 

N 

0113 

OEF 

•  ♦1 

a 

0114 

4MP 

DR  VR .  I 

0113 

» 

on* 

f 

ROUTINE 

TO  UA  IT  THOU  TINES  FOR  FLAG 

0117 

IMP 

NOP 

4119 

4  HP 

LIST 

0124 

LOA 

TWO 

0121 

STA 

IFLAG *  I 

0122 

CLC 

LDV 

0123 

CLP 

LDV 

4129 

J  SB 

9LIBX 

0123 

DEF 

♦  ♦l 

0120 

OEF 

•  ♦1 

4127 

4  HP 

DR VR  ,  I 

4  128 

* 

0129 

* 

0130 

*  VARIABLE 

STORAGE 

0131 

DC 

OCT 

60000 

0132 

K 1 

DEC 

312 

0133 

ZERO 

DEC 

0 

0134 

ONE 

DEC 

1 

0133 

TUO 

DEC 

2 

0  138 

THRE 

DEC 

3 

0137 

THOU 

DEC 

-32000 

0138 

CHT1 

NOP 

0139 

CHT2 

NOP 

0M4 

CHT3 

NOP 

0141 

PTRl 

NOP 

0142 

PTR2 

HOP 

0143 

VA  IT 

NOP 

0144 

END 
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APPENDIX  B  -  MUN  PCKXBAN  PROGRAM 


6L0V  T-00003  IS  ON  CO 40 0  41  USING  00094  BLKS  R«0000 
0001  FTN4  ,  L 

0002  C - C 

0043  C  C 

0004  PROGRAM  L0V!39l).  REV  830720'  TO  USE  THE  LDV 

0403  C  C 

0006  C  URITTEN  BY  STEVE  LUST1G  AND  CAROLYN  KAPLAN  C 

0407  C 

0008  C - FUNCTIONS  - - - - 

0009  C  1)  TO  TEST  THE  COINCIDENCE  FILTER  BOARD 

0010  C  2)  TO  OPERATE  THE  DISA  LASER  DOPPLER  ANEMOMETER 

0011  C  SYSTEM  HARDWARE  FROM  HP-2  IHX  SOFTWARE 


0012  C  3)  TO  COLLECT  DATA  FROM  DISA  LASER  DOPPLER  ANEMOMETER 
0013  C  BUFFER  INTERFACE 

0414  C  4)  TO  CREATE  OISK/TAPE  FILES  OF  THIS  DATA 

0013  C  3 )  TO  ASSEHBIE  DATA  t  WRITE  OUT  OH  LP  (LU6) 

001 t  C  4)  TO  COMPUTE  REYNOLD'S  STRESS  CORRELATIONS  ! t  ETC) 
0417  C  FROM  ASSEMBLED  DATA 

0018  C 


0019  C - MISCELLANEOUS - - - - 

0420  C  1)  THE  DRIVER  OPERATING  THE  LDV  IS  CURRENTLY  A 
0421  C  RELOCATABLE  VERSION  WHICH  IS  LOADED  WITH  THIS  PROCRAH. 
0022  C  THE  VARIABLES  PASSED  TO  THE  DRIVER! IDRVR  )  MUST  NOT  BE 
0423  C  CHANCED. 

0424  C  2)  THE  LAST  WORD  IN  ARRAY  SCW  MUST  BE  THE  'STATIC  CONTROL 
0023  C  WORD'  WHICH  IS  USED  TO  'TURN  ON'  DESIRED  CHANNELS.  PRE 
0026  C  CEEDING  CONTROL  WORDS  OPERATE  PERIPHERAL  BOARDS/DEV . S . 
0027  C  3)  THE  SORTING  ROUTINE  IH00E-3J  WILL  HOST  LIKELY  HAVE  TO 

0028  C  BE  MODIFIED  AS  MORE  HARDWARE  IS  ADDED  TO  THE  LDV  SYS 

0029  C  4)  RUHNING  THIS  PROGRAM  WITH  THE  FTN4.D  OPTION  WILL  ACTI-  C 
0430  C  VATE  PRINTING  SUBROUTINES  IH  TASK  3. 

0031  C . — . - . C 

0032  C 

0033  COMMON  LUO) 

0034  INTEGER  I  SC W< 1 0 ) , I D AT A< 3 1 2 > , I FL AC , T I TL E< 20  )  , IDCB (  636 > 

0033  INTEGER  T APELU , I T I H E< 3 ) , IBUFF! 2 0 ) , IDCB2( 1 44  )  , N AME( 3 > 

0436  INTEGER  NAM £2! 3>,PL0C,PARTCH, ST IME , ICNT /WAIT 

0037  REAL  FTI.FTW. FTC, TTI.DT. COUNT, F0RM1.FQRN2.RDATA 

0038  DATA  TITLE/40*2H  / 

0039  C 

0040  C - DETERMINE  THE  LU  OF  THE  DEVICE - - - 

0041  C 

0042  CALL  RHPAR(LU) 

0043  C 

0044  C - DETERMINE  TASK . . . . . . . 

0043  C 

0046  2  WRITE! LU, 3) 

0047  3  FORMAT! IX, ' ENTER:  1  TO  TEST  LDV  BY  TAKING  SAMPLES' 

0048  ♦/17X, ' (WRITE  TO  ANY  NON-OISC  LU)' 

0049  ♦711X»'2  TO  COLLECT  DATA  FROM  THE  LDV  INTERFACE' 

0030  ♦/ 1 7X . ' ! COLLEC T  FOR  A  GIVEN  AMOUNT  OF  TINE,' 

0031  4/17X, 'WRITE  RAW  DATA  TO  A  CREATED  DISC  FILE)' 

0032  ♦/ 1 1 X , ' 3  TO  SORT  AND  ASSEMBLE  COLLECTED  DATA' 

0433  +/17X, ' (READ  RAW  DATA  FROM  ANY  DISC  LU, ' 

0434  4/17X. ' ASSEMBLE  DATA,  WRITE  TO  DISC  FILE)' 

0033  ♦/ 1 l X , ' 4  TO  COMPUTE  TURBULENCE  STRESS  CORRELATION' 

0036  +/17X, ' VALUES  FROM  ASSEMBLED  DATA  FILES' 

0037  ♦/ 1 1 X , ' 3  TO  END  THIS  PROGRAM' //l  IX, ’ _'  ) 

0038  READ!LU. *  )  INODE 
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nnoonnoonononnoooooon 


0  95* 
0060 
0  06  t 
0062 
0063 
0064 
0063 
0066 
006? 
0068 
0069 
00?0 
0071 
0072 
0073 
0074 
0075 
0076 
0077 
0079 
0079 
0080 
0081 
0082 
0083 
0084 
0083 
0086 
0087 
0088 
0089 
0090 
0091 
0092 
0093 
0094 
0093 
0096 
0097 
0098 
0099 
0100 
0101 
0102 
0103 
0104 
0103 
0106 
0107 
0108 
0109 
0110 
0111 
0112 
0113 
0114 
0113 
0116 
Oil? 
0118 


CO  ro  i too . 200,300 , 400, 9999 >, INODE 

c 

100  CONTINUE 

C - TEST  THE  SYSTEM  AND  WRITE  TO  ANY  NON-DISC  LU - - --C 

C 

0  ASK  FOR  AMOUNT  OF  DATA  AND  OUTPUT  LU 
ISL-0 
IDL-512 
ICOUNT-O 
UR  I T  E( LU-  1 0  9  > 

109  FORMAT < '  ENTER  •  DATA  WORDS  TO  COL LECT Y <-3 1 2 )C 13 3 1  ) 

READY  LU- *  >IDL 
I F<  IDL  CT  312  UDL-312 
URITEY LU. 130) 

130  FORHATY'  ENTER  OUTPUT  LU<I2>s  -'  > 

READYLU, 131  )I OUT 

131  FORM ATY 12  ) 

C 

C  READ  STATIC  CONTROL  UORDS  AND  ZERO  BUFFER 

102  URITEY LU. 103) 

103  FORHATY'  ENTER  A  CONTROL  UORD  <000000  TO  STOP)  C  0  6  3 :  .') 

READ' LU, 104  )I£HTER 

104  FORMATY06) 

IFY I  ENTER  .  EQ .  08)  GO  TO  105 
I SL ■  ISL*  1 
ISCUY ISL >  *  I  ENTER 
GO  TO  102 

103  DO  110  1-1,512 

110  IDATAYD-0 
IFY  ISL  EQ  0  )  GO  TO  2 

C 

I  CODE-3 

CALL  DRVRY I  CODE,  IDATA, ISCU, ISL. ICNT. IFLAG) 

C 

C  WRITE  OUT  ANY  DATA  COLLECTED 

IFY ICNT. GT. 0)  UR  I TEY l OUT . 1 08  )  < ID A TA( I  ) , I • l ,  ICNT  ) 

108  F0RHATY63</8< IX, 06))) 

URITEY IOUT,  120  )  ICNT, IFLAG 
120  FORHATY'  ICNT  -  ',13/'  IFLAG  -  '.I2//> 

C 

C  COLLECTED  ENOUGH  ? 

ICOUHT-ICOUNT+ICNT 
IFY ICOUNT  .LT. IDL  )  GO  TO  103 
GO  TO  2 
C 

200  CONTINUE 

C - COLLECT  DATA  FROM  LDV  INTERFACE - - - 

C  SELECTIONS  OF  STATIC  CONTROL  UORDS  CAN  BE  ALTERED/STORED  HERE 
C  FOR  ROUTINE  DATA  COLLECTIONS.  THE  DISA  MANUALS  SHOULD  BE 
C  CONSULTED  BEFORE  MAKING  HARDWARE  CHANGES  ANY  HARDWARE  HODI- 
C  FI  CAT  ION  U ILL  PROBABLY  RESULT  IN  THE  NEED  TO  MODIFY  CONTROL 
C  UORDS.  THIS  IS  A  GOOD  PLACE  TO  STORE  CONTROL  WORDSi 
ISCUY  1  >-0  1001  IB 
ISCUY2  )" 0  22  34  2  B 
ISCUY3  >-0000348 
ISL-3 
IDL-312 
UAIT-5 

C  1  ISCUY  1  )  TURNS  OH  THE  COINCIDENCE  FILTER  BOARD  AND  CHANNELS  C 
C  ONE  AND  TWO.  C 


4 


'  i 

•  1 


•  4 


I 
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Oils  C  2.  ISCW( 2 )  SETS  A  125  HICROSEC  TIME  INTERVAL  FOR  A  0.1 

0120  C  HICROSEC  TIHE  BASE  FOR  THE  COINCIDENCE  FILTER  BOARD 

0121  C  3  ISCW<3>  SAMPLES  BOTH  PART  CHAHHELS  OF  CHAHHELS  1.2  AHD  4. 
0122  C  CORRESPONDING  TO  P1.P2  AND  P4  ON  THE  INPUT  MULTIPLEXER 
0123  C  BOARD. 

0124  C  4.  THE  DRIVER  ASSURES  THAT  THE  ADDRESS  TO  OROP  REQUEST  B  IS 

0123  C  SET  AS  SB  <  THE  SWITCH  ON  THE  57GXXX  BORD  IS  OH  6> 

0126  C  TAKE  OATA  CALL  I HG  USER  WRITTEN  DRIVER  AND  PASS  THE  FOLLOW  IHG 
0127  C  VARIABLES: 

0128  C  ICODE  :  -1  PERFORM  INITIATION  SECTION  OF  DRIVER  ONLY 

0129  C  (CLEAR  BUFFER  MEMORV  AND  SEND  CONTROL  WORDS) 

0130  C  -2  PERFORM  CONTINUATION  SECTION  ONLY 

0131  C  (COLLECT  ANV  DATA  IN  BUFFER  MEMORY) 

0132  C  *3  00  BOTH  SECTIONS  WITHOUT  RE TOR HI HG  TO  THIS 

0133  C  PROGRAM  IN  BETWEEN 

0134  C  ID ATA( 5 1 2  )  :  INTEGER  ARRAY  OF  DATA  TRANSFER  WORDS  COLLECTED 

0133  C  IDL  :  INTEGER  HUMBER  OF  DATA  SAMPLES  TO  COLLECT( »<5 12 

0136  C  ISCW(IO)  i  INTEGER  BUFFER  OF  CONTROL  WORDS  TO  SEND  TO  THE 

0137  C  57G20  CONTROL  SYSTEM. THE  LAST  WORD  SHOULD  BE 

0138  C  THE  STATIC  CONTROL  WORD  (SCW). 

0139  C  ISL  :  INTECER  HUMBER  OF  CONTROL  WORDS  TO  BE  SENT. 

0140  C  ICNT  :  •  OF  DATA  WORDS  COLLECTED  BY  THE  DRIVER 

0141  C  (INITIATED  AT  THE  BEGINNING  OF  EVERY  CALL) 

0142  C  IFLAC  -1  MEANS  THAT  THE  DRIVER  COLLECTED  DATA  NORMALLY 

0143  C  -2  THE  DRIVER  HAD  TO  WAIT  TOO  LONG  FOR  THE  FLAG 

0144  C  FLIP-FLOP  TO  SET  WHILE  RTE  WAS  SUSPENDED. 

0143  C  NOTE  THAT  IDATA  BUFFER  IS  NOT  CODED  ON  DISC  IN  ASCII  CHARCTR 

0146  C  SET.  THUS  IT  IS  IMPOSSIBLE  TO  SEE  THIS  DATA  W/'OUT  SOFTWARE. 

0147  C  IN  ADDITION.  THE  DATA  BUFFER  IS  HOT  2ER0ED  BETWEEN  DR VR  CALLS 
0148  C  TO  SAVE  TIHE  THIS  PROGRAM  KNOWS  HOW  MANY  HEW  WORDS  HAVE  BEEN 
0149  C  COLLECTED  BY  THE  COUNTER  ICNT.  THUS  ONLY  THE  FIRST  ICNT  WORDS 

0130  C  IN  BUFFER  IDATA  HAVE  BEEN  COLLECTED  AFTER  ANY  DR VR  CALL! 

0131  C - - - - - . 

0132  COUNT-O.O 

0133  DATA  T I TL E/40* 2H  / 

0134  C 

0133  C  GET  INFORMATION  FROM  USER 
0136  WRITE(LU.2tl) 

0137  211  F ORH AT ( 3X . '  ENTER  CARTRIDGE  8.  OUTPUT  FILE  NAME! 12. A6 I  .') 

0138  READ< LU . 202  )  ICR. NAME 

0139  202  FORHAT( 12. IX. 3A2  ) 

0160  WRI TE( LU . 204 ) 

0161  204  F0RHAT(3X. '  ENTER  APPROXIMATE  TINE  IN' 

0162  +/6X ,  ’  SECONDS  TO  COLLECT  DATA(UO):  .'  ) 

0163  READ( LU. 208  )  STINE 

0164  208  FORHAT(IIO) 

0163  WRITE( LU. 209) 

0166  209  F0RHAT(3X, '  ENTER  A  TITLEIA361') 

0167  READ(LU. 210  )( TITLE( I), I *2, 20) 

0168  210  FORMAT ( 20 A2  ) 

0169  C 

0170  C  CREATE  AND  OPEN  A  TYPE  3  SEQ  ACCESS  FILE 

0171  ISIZE--1 

0172  IOPTH-0 

0173  I T YP  *3 

0174  I  SEC U* 0 

0173  IDC8S-640 

0176  CALL  CREAT( IDCB. IERR, NAME. ISIZE.ITYP.ISECU, ICR. IDCBS) 

0177  I  F(  I  ERR .  LT  .  >)  >  CALL  ERR(  1  .  I  EPR  .  I  ) 

0178  CALL  OPEHdDCB,  IERR,  NAME.  IOPTN,  ISECU,  ICR.  IDCBS) 


oooor/oonooonoooowoooonnnoooo  non 


0171 

0194 

0191 

0192 

0193 

0194 

0193 

0186 

0197 

0199 

0199 

0190 

0191 

0192 

0193 

0194 

0193 

0196 

0197 

0199 

0199 

0240 

0201 

0202 

0203 

0204 

0203 

0206 

0207 

0209 

0209 

0210 

0211 

0212 

4213 

0214 

0213 

0216 

0217 

0219 

0219 

0220 

0221 

0222 

0223 

0224 

0223 

0226 

0227 

0229 

0229 

0230 

0231 

0232 

0233 

0234 

0233 

0236 

0237 

0239 


i 


< 


IF!  I  ERR  L 7 . 0  >  CALL  ERR! 2  < I  ERR  <  2 ) 

C 

C  WRITE  OUT  A  TITLE 
NU0RDS-20 

CALL  URITF< IDC8, IERR. TITLE- NUQRDS  )  { 

IF< I  ERR. LT  .  0  )  CALL  ERR<  3, IERR,  3  ) 

C 

C  GET  STARTIHC  TINE 

CALL  £XEC(11<  I TI HE  .  I  YEAR  ) 

FTI-! ITIHE<  4>*60*ITINE<  3  ))"60*ITINE! 2) 

C 

C  COLLECT  AND  WRITE  OUT  A  BUFFER  OF  DATA  -  ' 

203  I  CODE* 1 

CALL  DRVRUCODE. IDATA.IDL. ISCW, ISL . I CN T . I FL AG > 

CALL  EXECOl.  ITINE,  IYEAR) 

TTW" !  ITIHE(  4)*60+ITIME<  3  ))»60+ITIHE!  2  > 

203  CALL  EXEC! 1 1, ITINE. IYEAR  ) 

NTM*WAIT*< ITINE! 4)*60+ITINE! 3D  +  ITIHE!  2  ) 

I F! NTW . LT  .  TTW  )  GO  TO  203 

ICODE-2  < 

CALL  DRV R< ICODE<IDATA<IDL<  ISCU<ISL<ICNT<  I FLAC> 

C 

C  WRITE  BUFFER  TO  TERMINAL,  IF  DESIRED 
0  IF<  I CNT . G  T . 0  >  UR  ITE<  LU .  108  )<  I  DATA<  I), 1*1,  ICNT> 

D  WRITE! LU.  120)  ICNT <  IFLAG 

C 

C  WRITE  OUT  DATA  TO  DISC  FILE  , 

IF< ICNT. GT. 0>  CALL  UR  I TF< I DCB < I  ERR < I  DAT A ,  ICNT > 

I F< I  ERR  LT  0)  CALL  ERR< 3 , 1  ERR < 3  ) 

COUNT* COUNT ♦ICNT 
C 

C  CHECK  TIME 

CALL  EXECdl.  ITINE.  IYEAR) 

TTI ■< ITIHE<  4)* 60  + 1  TINE!  3  )  )*60  +  I TI HE< 2 > 

DT-TTI-FTI 

IF!DT.LT. STINE)  CO  TO  203 
RATE "COUNT/DT 
C 

C  WRITE  OUT  •  SAMPLES  AND  CLOSE  THE  DISC  FILE 
URITEILU, 207)  COUNT . RATE , NANE < I CR 
207  FORH AT! I  OX . '  SAMPLING  TINE  OVER!' 

♦/10X.G13 . 3. '  DATA  WORDS  COLLECTED ' 77X. ' AT  ',G13.3< 

♦  '  WORDS/SEC ' /7X ,  'ON  '  <  3A2.  '  s  '■  ' .  12/7  ) 

CALL  LOCFUDCB<IERR<IREC<IRB<  IOFF,  JS  EC) 

IF!  IERR. LT. 0)  CALL  ERR! 3 . 1  ERR < 6 > 

ITRUH«JSEC/2-IRB-l 

CALL  CLOSE!  IDCB.  IERR.  ITRUN  ) 

IF!  IERR. LT . 0)  CALL  ERR! 4 , I  ERR < 7  ) 

C 

C  GO  BACK  TO  COMMAND  SECTION 
CO  TO  2 

C 

C - SORT  6  REDUCE  RAW  DATA  FILES - - - 

300  CONTINUE 
C 

C  PREPARE  INPUT  OISC  FILE  WITH  RAW  DATA 

I OPTN" 0  _  < 

ITYP-3 
I S I 2E»- 1 
ISECU-0 


I 
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1239 
■  :->» 
0  2  4  1 
0242 
0243 
0244 
0243 
0244 
0247 
0248 
0248 
0250 
0231 
0252 
0233 
0234 
0233 
0234 
0237 
0238 
0233 
0240 
0241 
0242 
0263 
0264 
0263 
0264 
0267 
0268 
0269 
0270 
0271 
0272 
0273 
0274 
0273 
0276 
0277 
0278 
0279 
0280 
0281 
0282 
0283 
0284 
0283 
0286 
0287 
0288 
0289 
0290 
0291 
0292 
0293 
0294 
0293 
0296 
0297 
0298 


I  DC  8  S*64  0 
IDCBS2-128 
WRITE'.  LU  -  3«4> 

3<M  F ORMAT  <  3X  ,  '  ENTER  CARTRIDGE  *  ,  INPUT  FILE  NAME  CI2.A61  .') 
READ! LU, 202  )  ICR, NAME 

CALL  OPEN  <  I  DCS  .  IERR  -NAME  .  I  OP  TN .  ISECU ,  I  C« ,  IOCBS  ) 

IF<  I  ERR . L  T . 0  >  CALL  ER  R<  2  ,  I  ER  R  ,  8  ) 

C 

C  CREATE  AND  OPEN  OUTPUT  FILE 
URITE<  LU, 21  1  > 

READILU. 202  )I0UT,NAHE2 

CALL  CREAT< IDC82, IERR,NA«E2. I  SIZE, IT  VP, ISECU, I  OUT, IDCBS2  > 
IFMERR.  LT.  0>CALL  ERR( 1 , IERR , 9 > 

CALL  OP£N< I DCB2 , IERR, NAHE2- I OPTN .  I SE CU , 10 UT ,  IDCBS 2  ) 
IF<IERR.LT  0>  CALL  ERR<  2 , I  ERR , 1 0  > 

c 

C  WRITE  ASSEMBLED  DATA  TO  A  NON-DISC  LU ,  IF  DESIRED 
0  «RITE< LU, 305) 

0303  FORM AT<  SX . '  ENTER  A  NOH-OISC  LU  TO  SEE  ASSEMBLED  DATA:  .') 
D  READ<  LU, 131  )  I LU 

D  URITE( LU, 30b)  ILU 

0306  FORMAT  <  3X , '  DO  YOU  HEED  TO  CONTROL  LU  ',t2.'  7CY7H1  _') 

D  READ( LU, 307  )  IANS 

0307  FORMAT ( A  1  ) 

0  IFMANS.EO  1HY)  CALL  TAPE<ILU) 

C 

C  READ  IN  TITLE  AND  WRITE  TO  OUTPUT  FILE 
NW0RDS-20 

CALL  READFI IDCB, IERR,  IB UFF,N WORDS .  ILEN  > 

IFMERR.LT. 0)  CALL  ERR<  4  ,  I  ER  R ,  1 1  > 

CALL  WRITFC IDCB2, IERR- I  BUFF. NWORDS  ) 

IFMERR.LT.  0>  CALL  ERR<  3 , 1  ERR  ,  l  2  > 

NCHRS-40 

CALL  CODE(NCHRS) 

READ< I  BUFF, 2tO  X  TITLE<  t  >, t«t . 20  ) 

C 

C  INITIATE  STATISTICAL  QUANTITIES.  IPASS-1 
IPASS-1 

CALL  STAT( I  PASS, ICNTR.RO AT A,  I  TYPE  ) 

C 

C  CALL  FIRST  PRINTING  ROUTINE,  IF  DESIRED 
D  CALL  PRIN1< TITLE. ILU  ) 

C 

C  RE  AO  IN  A  RECORD  OF  DATA  OF  UP  TO  312  WORDS 
HCHRS-22 
HWRDS-11 
NWORDS-312 
I  PASS*  2 
309  CONTINUE 

DO  308  1-1.512 
303  l DAT A<  I  )-0 

CALL  READF< IDC 8 . IER R .  ID A T A , N WORDS ,  IL EH  ) 

IFMERR.LT. 0>  CALL  ERR<  6  ,  I  ERR  ,  14  > 

IF< ILEN. LT. 0)  CO  TO  311 

C 

C  PROCESS  THE  DATA  ON  A  RECORD 
ISKIP-0 

DO  310  1*1, ILEN 
C 

C  ASSEMBLE  THE  DTW  INFORMATION 
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0299  CALL  ASM8L<  I0ATA<  I  ) , PARTCH, PLOC , FORM t - F0RM2,  IBO , I  TO  ) 

0  3  0  0  : 

030  1  •:  SORT  THE  INFORMATION  IN  THE  DTK 

0  30  2  CALL  SORT< IOATA<  I  ), PARTCH. PLOC-  I  TYPE. ICNTR.F0RM1, F0RM2.RDATA) 

0303  I F< I  TYPE . EG  6 >  GO  TO  310 

o  3  4  : 

0305  :  PRINT/WRITE  THE  DATA.  IF  DESIRED 

0304  0  CALL  PRIN2MCNTR,  ITYPE.RDATA,  ISO,  ITO,  ILU) 

0  30  7  C 

03oS  C  DO  STATISTICAL  SUMMATIONS,  IPASS«2 

0309  CALL  STAT< IPASS, ICNTR, ROATA, ITYPE) 

0  310  C 

0311  C  CROUP  VELOCITY  i  INTRVL  DATA  (PER  BURST  )  l  KRITE  TO  OUTPUT  FILE 

0312  CALL  GATHR<  I SK IP  ,  ICNTR, ITYPE, RDATA.IDCB2) 

0313  310  CONTINUE 

0  3  t  4  : 

0315  C  READ  THE  NEXT  RECORD  IK  THE  FILE 
0316  GO  TO  309 

0317  C 

0318  C  DO  FINAL  STATISTICAL  COMPUTATIONS  ON  A  BUFFER  OF  DATA,  IPASS=3 
0313  311  IPA5S»3 

0320  CALL  STATC IPASS,  ICNTR, RDATA. ITYPE  ) 

0321  C 

0322  C  CLOSE  THE  OUTPUT  FILE 

0323  CALL  LOCF< IDCB2,  IERR, IREC, IRB, IOFF, JSEC  ) 

0324  IFMERR.LT. 0>  CALL  ERR< 3 , l ERR , l 6 > 

0325  ITRUH«JSEC/'2-IRB-l 

0326  CALL  CLOSE< IDCB2, IERR, ITRUN) 

0327  IFMERR.LT. 9)  CALL  ER R<  4  ,  I  ERR ,  1 7  > 

0328  C 

0329  C  PLACE  EOF  ON  ILU  IF  IT  IS  A  TAPE 

0330  ICHHDMLU  +  100B 

0331  CALL  EXEC(3,ICNUD>0) 

0  33  2  C 

0333  C  RETURN  TO  COMMAND  SECTION 
0334  GO  TO  2 


0335  C . - . - . - . . . . C 

0336  400  CONTINUE 

0337  C - COMPUTE  REYNOLDS  STRESSES . - . - . . C 

0338  9999  URITE<  LU. 9998  ) 

0339  9998  F ORM AT ( I  0 X . ' *  *  END  OF  LDV«*'> 

0340  END 

0341  c**************«************************************************C 


0342  SUBROUTINE  SOR T<  I D A T A , P ARTCH , PL OC .  I T YP E , I CNT R , FORM  1  , 

0343  +  F0RM2, RDATA  ),  TO  DECIPHER  DTK'S 

0344  COMMON  LU<5) 

0345  INTEGER  P 1 2 34 < 4 ) , D I G I 0< 4  ) . HOD E< 4 > , PARTCH , PLOC 

0346  REAL  F SH I  FT ( 4 ) , L AMB D A< 4 ) , THE T A< 4 > , FORM  I , F ORM 2 , TB A SE , RD AT  A 


0347  C . . . C 

0348  C* ♦**  *  NOTICE  TO  USER  MAKING  MODIFICATIONS  **•*•  C 
0349  C  INDICATE  DEVICES  CONNECTED  TO  THE  57G120  INPUT  MULTIPLEXER  s  C 
0350  C  THE  USER  TELLS  THIS  PROGRAM  HOW  THE  57G20  CABLES  ARE  CON-  C 
0351  C  HECTED  BY  SETTING  THE  MEMBERS  OF  THE  FOLLOKING  ARRAYS:  C 
0352  C  P1234(X)-V  WHERE  X«l,2,3,  OR  4  C ORRE SPO HD  I  MG  TO  PINS  C 
0353  C  P l . P2 , P  3 , P  4  RESPECTIVELY.  C 
0354  C  Y-1,2  FOR  COUNTER  WHICH  WILL  BE  REF-  C 
0355  C  EREHCED  AS  1.2  C 
0356  C  Y«9  FOR  THE  COINCIDENCE  FILTER  BOARD  C 
0357  C  Y»0  IF  NOTHING  IS  CONNECTED  AT  ALL  C 
0353  C  OR  FOR  THE  COINCIDENCE  FILTER  BOARD  C 


0339 
0360 
0361 
0362 
0363 
0364 
0365 
4366 
4  36  7 
0369 
0369 
0370 
0371 
0372 
0373 
0374 
0375 
0376 
0377 
0378 
0379 
0380 
0381 
0382 
0383 
0384 
0385 
0336 
0387 
0388 
0389 
0390 
0391 
0392 
0393 
0394 
0395 
0396 
0397 
0  398 
0  j9  9 
0400 
040  1 
0402 
0403 
«  4  •,  4 
0  405 
0406 
0407 
0408 
0409 
0410 
0411 
0412 
0413 
0414 
0413 
0416 
0417 
0418 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

5 

c 

c 

c 

c 

c 


0 1 G I  0 ( X  )■  A  WHERE  X  IS  THE  SAME  A  S  ABOVE  C 

A»1  IF  THE  DIGITAL  I  70  SLOT  41  CORRES  C 
TO  PIN  X  OR  A *2  IF  DIGITAL  I/O  SLOT  C 
*2  APPLIES  C 

A-0  IF  A  COUNTER  IS  HOT  CONNECTED  C 

FSHIFT(X>*8  WHERE  X  IS  THE  SAME  AS  ABOVE  C 

B  IS  THE  FREQUENCY  SHIFT  AS  SHOWN  OH  C 
THE  55N14  FREQ  GENERATOR,  IE  6  IS  C 
POSITIVE  IF  THE  >49MH2  SIDE  IS  USED  C 


AND  NEGATIVE  IF  :4oMH2  B  IS  ONLY  !K£C 
VALUE  SET  ON  THE  55H14  PANEL  C 

«ODE(X)*C  WHERE  X  IS  THE  SAME  AS  ABOVE  C 

C=1  IF  DATA  IS  COLLECTED  IN  THE  FIXED  C 

MODE  ON  THE  COUNTER  PROCESSOR,  C 

C-2  IF  THE  COMBINED  MODE  IS  USED.  C 

C*  3  IF  THE  VARIABLE  FRINGE  MODE.  AND  C 

C-4  IF  THE  TRANSIT  MODE  IS  USED  C 


C«9  IF  HO  COUNTER  IS  OH  X  .  C 

TBASE  IS  THE  FLOATING  POINT  CLOCK  TIME  BASE  C 

OH  JUMPER  S  2  OF  THE  CFB  IT  MAY  BE  0  1  C 
1 ,  OR  10  USEC  C 

LAMDA(X)-D  WHERE  X  IS  THE  SAME  AS  ABOVE  C 

D  IS  THE  WAVELE  HGTH  OF  LASER  EUR  USED  C 
OH  CHANNEL  X.  C 

THETA(X)«E  WHERE  X  IS  THE  SAME  AS  ABOVE  C 

E  IS  THE  ANGLE  BETU/  THE  CORRES  BEAMS  C 
ASSOC  W/  THIS  COUNTER  (FOR  COAXIAL  C 
BEAMS  THIS  IS  ALUAYS  CONSTANT)  C 


DATA  P  123  4/ 2 ,  1.0,9/ 

DATA  D  IG  I  0/2, 2,0,0/ 

DATA  FSHIFT/-10  OE  +  03,-10  OE ♦ 0 3 , 0  0 , 0 .  0  / 
DATA  MODE / l ,1,9,9/ 

DATA  T8ASE/0 . 1 E-06/ 

DATA  LAMBDA/  4 8 8  0 E - 0 9 ,  5  1 4  5 E -0 9 , 0 . 0 , 0  0/ 
DATA  THETA/O. 019330253, 0 .019701965 ,0 .0, 0  0/ 


ITYPE  IDENTIFIES  THE  TYPE  OF  INFORMATION  THAT  A  DATIJM  GIVES:  C 


I T  Y  PE  *  1  FOR  DOPPLER  FREQUENCY  DATA  C 

2  FOR  SAMPLE  INTERVAL  DATA,  C 

3  FOP  HUMBER  OF  FRINGES,  C 

l  FUR  BURST  TIME  DATA,  C 

5  FOR  VELOCITY  DATA,  AND 

6  I  DATA--') 


C  IDEHTIFY  THE  TYPE  OF  DATA 

ITYPE»DIGIO(PLOC  )*PARTCH 
I F<  ITYPE  EQ  2  AND . PAPTCH  E Q  IB>  ITYPE-3 
C  SAMPLE  INTERVAL  HhS  ONLY  ONE  PART  CHAHHEL  OF  INFO: 

I F<  ? 1234c  PLQC  i  EQ  9  AND . PART  CH  EQ  IB)  ITYPE«2 
I  F<  P  1234  c  PLOC  i  EQ  9  AND  PARTCH  EC!  2B  >  ITYPE  =  6 
I  F-  IDuTA  EQ  0  >  I  TYPE  =  6 
C 

C  CONVERT  DIGITAL  OUTPUT  TU  DOPPLER  FREQ,  BUR  ST  ,’T  R  AN  S I  T  TIME, 

C  HUMBER  OF  FRIHGES,  OR  SAMPLE  1N1ERVAL  TIME  THIS  IS  DONE  AS 
C  PRESCRIBED  OH  P27  OF  DISA  COUNTER  INSTRUCTION  MANUAL.  ALSO 
C  SUBTRACT  FREQUENCY  SHIFTING  AND  CONVERT  TO  FREQUENCY  TO  VELOCITY 
IF'.  MODE!  PLOC  >  GT.2>  GO  Tr)  303 
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>n;*  ;  f  1  itype  eo  h  p  dh  ta=fgrh  i  *1 4  96+fshift«:  ploo 

••420  IF'.  ITYPE  E  Q  4)  R  DA  Th  «  1  6  E  +  09 /F  0  RM  1  ♦  F  SH  I F  T  (  P  LOC  ) 

.  -4 ;  •  IF<  ITYPE  £ Q  4  )  I  TYPE"  1 

j  4 .  ”5'  T  c'  MODE!  PLOC)  LE  2)  GO  TO  354 

0423  'C.  !TVP£  EQ  l  )  RDATA-0  53  4  82  FO  R  K  1  *F  S  H I  FT  <  PL  OC  > 

0424  l  F •  ITYPE  EQ  4)  R DR TA  =  0  5E - 09 *F0 RM 1  +  FSH IF T ( PLOC  ) 

,i42  5  I  F«  ITYPE  E8  4  >  ITYPE-  i 

0426  C  CONVERT  DOPPLER  FREQUENCY  TO  VELOCITY 
0427  354  IF.  ITYPE  NE  l  )  GO  TO  355 

0428  C'D>  t  F.  t  £8 .  1  )  GO  TO  355 

0  424  FACTOR-LAMBDAt  PLOC  »/<  2. *SIH<  THETA< PLOC  1/2  .  > > 

>43)  RDATft-  ROhTN*FhCTOR 

0431  ITYPE-5 

0432  355  I  Ft  ITYPE  EO  21  RDATA*  F0RM1*TBASE 

0433  IF'.  ITYPE  EQ  3)  R  D  A  T  A  3  FORM  2 

0434  ; 

0435  C  IDENTIFY  WHICH  COUNTER 

0436  ICNTR=P1234<PlOC 1 

0437  C 

0438  C  LEAVE 

0438  RETURN 

0440  £  HD 

0442  SUBROUTINE  GAT  HR < I SK  I  P .  I  CHTR  .  IT YPE , RDAT A , IDC B2  )  .  TO  CROUP  DATA  BY 

0443  ♦  8UR3TS 

0444  COMMON  L U ( 5  ) 

0445  DIMENSION  I BU F F< 33  )  , I T YP ( 1 0  )  •  ICNT (  1 0  )  , R D AT<  1  0 > 

0446  I ADDRS-I CNTR* I  TYPE 

0  44  7  C 

0448  C  IF  THIS  IS  THE  FIRST  CALL.  SAVE  THE  FIRST  ADDRESS 
0449  I  Ft  I SK  IP . EQ . 0  ICO  TO  506 

0  45  0  : 

0451  C  IF  IT  IS  NOT  THE  FIRST  CALL,  CHECK  TO  SEE  IF  ADDRESS  EQUALS 
0452  C  FIRST  ADDRESS 

0453  IF<  I ADDRS  NE.  I  AD  )G0  TO  522 

0454  C 

0455  C  IF  ADDRESS  EQUALS  FIRST  ADDRESS ,  ZERO  THE  ARRAY.  AND  RESET  IPTR 

0456  52  1  I P T R  *0 

0457  DO  508  1-1.3 

0458  ITYPCI1-0 

0459  ICHTiD-0 

0460  RDAT <  11-0  0 

0461  Sod  CONTINUE 

0  46  2  C 

0463  C  IF  DATA  IS  VELOCITY  OR  SAMPLE  INTERVAL  DATA  WRITE  IN  ARRAY 

0464  522  IF«  ITYPE  HE  2  AND  ITYPE  NE  5)  RETURN 

0465  I PTR-I PTR  ♦  1 

0466  I SK I P  *  1 

0467  I DUM  *  I CN  T  R 

0468  IF'  ITYPE  .  EQ  .2  )IDUM=>3 

0469  ITYP( IDUM  )»ITVPE 

0470  I  CNT< I DUM  )- ICNTR 

047  1  RD  A  T  < IDUM  1- RD A  T  A 

0  4’ 2  C 

0473  C  IF  ARRAY  IS  FULL.  WRITE  ARRAY  TO  BUFFER 
0474  IF.  IPTR  ME. 3 >  RETURN 

0475  CALL  CODE'.  N0HR3  < 

0  476  WRI TE<  IBUFF ,5ol  )< ICNT< IPTR  ).  ITYP< I  P  T  R  >  > RDATt  IPTR  ).  IPTR*!  .3  > 

0477  50  1  F0RMAT<3< IX. 12, IX, 12.  IX, G15  8  )> 

0478  C 


0479  C  WRITE  DA'A  IN  BUFFER  TO  OUTPUT  FILE  &  RETURN 
0480  CALL  WRI  TF<  IDCB2 , I  ERR  , I  BUFF , NWRDS  1 

0481  I  F<  I  ERR. LT. 0  )C ALL  ERR ( 5 ,  IE RR ,  1 5 > 

0482  RETURN 

0483  C 

0484  C  IF  FIRST  CALL.  SAVE  FIRST  ADDRESS  IN  VARIABLE  I  AD .  AND 

0483  C  THEN  2ER0  THE  BUFFER 

0488  306  IAD*IADRS 

0487  NCHRS-66 

0488  N UR D  S* 33 

0489  GO  TO  521 

0490  END 

0492  SUBROUTINE  STAT< IPASS . I CNTR, DATA,  I  TYPE >  .  TO  DO  STATISTICS 

0493  COM  NON  LU<5  ) 

0494  REAL  HAG . ANGL . SUHt , SUM2 . AVG1 . AVG2 

0493  INTEGER  N  1 , N2 

0496  GO  TO  <  100, 200.300  ),  IPASS 

0497  C 

0498  C  INITIATE  QUANTITIES 

0499  100  CONTINUE 

0500  N  1*0 

0501  N2*0 

0302  SUM  1 *0 

0503  SUH2*0 

0504  RETURN 

0505  C 

0506  C  DO  SUMMATIONS,  ETC 

0507  200  CONTINUE 

0508  IFMTYPE  HE. 5)  RETURN 

0509  IF<  ICNTR  .  HE  .  1  J  GO  TO  201 

0510  N  1»N  1*1 

0511  SUMi*SUMl+DATA 

0312  RETURN 

0513  201  N2*N2*  1 

0514  SUM2 "SUN2  +0 AT  A 

0513  RETURN 

0516  C 

0317  C  FINISH  CALCULATIONS 
0318  300  CONTINUE 

0519  A VG l *SUH 1  /N  l 

0520  A  VG  2  »SUM2 /N2 

0  321  MAG* SORT!  AVG1*AVG1*AVG2*AVG2  > 

0322  AHGL«ATAH<AVG2.'AVGM 

0523  ANGL»ANGL*130,'3  141592654 

0324  URITE<6.301)  A VG 1 . A VG 2 , M AG , AHCL 

0323  URITE< LU, 301  >  AVG 1 , AVG2 , HAG . A HGL 

0326  301  FORMAT!//'  UNWEIGHTED  MEAN  VELOCITIES:' 

0327  ♦  /'  <  CNTR  1>  '  ,015.8/'  <  CNTR  2>  '.G15.8 

0328  ♦//'  RESULTANT  VECTOR:  '/'  MAGNITUDE  *  '.G15  8 

0529  ♦/'  INCLUDED  ANGLE  «  '.G13  8//) 

0530  RETURN 

0331  END 

0332  C* ***••**•***••**•*****•*•*•****•*» ******* . . 

0333  SUBROUTINE  ASHBL ( D T W , PAR TCH , P LOC , F ORH 1 , F 0RM2 , 1  BO ,  IT 0  )  ,  ASSEMBLE 

0334  *  DATA  FROM  DTW'S 

0333  C  THE  FORMAT  OF  P.D,  I.  T  TYPE  DATA: 

0337  C  *  OTC  ♦  80  CCl  CC2  ♦  M7  E3  E2  ♦  El  EO  M6  ♦  M3  M4  M3  ♦  M2  Ml  MO  ♦ 
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>539 
0540 
0341 
0  54  2 
0543 
0344 
0543 
0346 
0547 
0548 
0549 
0550 
0351 
0352 
0353 
0554 
0555 
0356 
0557 
0558 
0559 
0560 
0561 
0562 
0563 
0564 
0565 
0566 
0567 
0568 
0569 
0570 
0571 
0572 
0573 
0574 
0575 
0576 
0577 
0578 
0579 
0530 
0531 
0532 
4533 
0584 
0535 
0536 
0587 
0588 
0589 
0590 
0591 
0592 
0593 
0594 
0595 
0596 
0597 
0598 


F7  F6  ♦  F5  F 4  F3  +  F2  FI  F0  ♦ 


I  BO i  ITO. [MAN. I  EXP 


C  *  *  *  •« 


:  o 
i  l 


TH  E  FORHAT  OF  FRINGE  DATA 
♦  OTC  ♦  80  CCl  CC2  ♦ 

f  ♦ 

COMMON  L U ( 5  ) 

INTEGER  DTU  .PARTCH, PLQC 
REAL  FOR  Ml. F0RN2 

PARTCH a  I AND(DTW,  1000006  )/ l 000008+ 18 
PLOC*IAND<DTU. 03  00  0 08  )710000B+1B 
I  80 alANDTDTH. 4000081. ’40  0008 
I  TO  a  IAHO>  DT  A  ,  4  •>■',  OD  >740008 
IMMHainNiu  D TO  177B  >+IT0*2008 
I  EXP  »  I  AND ( D TO • 36008  >7  2008 
F0RH1»IMAH*!2  **IEXP> 

F ORM 2=  IAN D!OTW, 37700)7108 
RETURN 
END 

**•**•*•«•*«*** 

SUBROUTINE  PR  I N l ! T I TLE,  ILU), 

COMMON  LU ( 5  ) 

INTEGER  I  SCW!  10>,TITLE<20).STARS<40> 

DATA  STARS 7 2H  39*2H**7 
URI TE<  ILU  .  1  1  >  (TITLE!  I),  1*1 
FORMAT-;  1  HI.  20 A27  IX.  COUNTER 

♦  6X .  'PART  CHANNEL' .  1  SX  > 'DATA  S ENT ' , 5X , ' 0 VFLOW  8 1 T '  ,  3X . ' T I  ME  OUT' 

♦  740A2  > 

RETURN 

END 

SUBROUTINE  PR  I H2 < CN TR , I T YP E . ROA TA .  IB  0 .  I  TO . I L U > .  SECOND  PRINT  ROUT; 

♦  I  NE 

COMMON  L  U ( 5  )  j 

INTEGER  CNTR, PARTCH. IBO. ITO, ALPHA! 15,5 > 

REAL  RDATA  > 

DATA  ALPHA72HD0. 2HPP. 2HLE, 2HR  , 2HFR, 2HEQ. 2HUE. 2HNC. 2HY  .6*2H  , 

♦  2HSA,2HMP,2HLE,2H  I  , 2HNT , 2HER , 2HVA , 2HL  ,7*2H  , 

+2H#  .2H0F.2H  F , 2HRI , 2HNG , 2HES , 9 *2H  ,  ; 

♦2HBU.2HRS.2HT  , 2HT I , 2HME . I 0*2H  .  i 

♦2HVE,2HL0,2HCI,2HTY,11*2H  7 

C  DO  NOT  PRINT  THE  SAME  DATUM  OVER  fc  OVER  ! 

C  IF! ROATA. EQ .PDATA>  RETURN 

C  P  DAT  A=»RO  A  T  A 

WRITE!  ILU,  1 0)  CNTR, ! ALPHA!  I .  I  TYPE  )  , I » l ,  15 > . ROATA , I B 0 . I  TO 
10  F  OP  M  AT  <  5  M , I  1 . 8  X ■  1 5  A  2 ■ G 1 5  .  4 , 9X , I  1 , 8  X ,  1 1  ) 

RETURN 

END 

SUBROUTINE  ERR! I  PASS  IERR ■  ILOC > .  TO  FIND  DISC  i  SOFTWARE  ERRORS 


FIRST  PRINTING  ROUTINE 


20  > ,!  STARS! I > , I ■  1 , 40  > 

*'  . 

SENT' ,  5X, ' OVFLOW  BIT' ,3X. 


COMMON  L U ! 5  ) 

IF! IPASS. EQ  1 > 
IF!  IPASS. EQ  2) 
IF*.  IPASS  .  EQ .3  > 
IF!  IPASS . EQ  .4  > 
IF!  IPASS  . EQ .5  > 
IF-.  IPASS.  EQ  6  > 
I  F>.  £  ERR  GE  .  0  > 
STOP 


URI TE<  LU , 
WRITE! LU, 
WRITE! LU , 
WRITE! LU, 
WRITE! LU, 
WRITE! LU, 
RETURH 


1 >IERR» 

2  HERR, 

3  > IERR , 
4 ) IERR , 
5 i IERR, 
6 ) IERR , 


ILOC 

ILOC 

ILOC 

ILOC 

ILOC 

ILOC 
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0999  l  FORHATc  '  <  TCREAT  HE  S  S A  GE  >  > l ER  R  *  AT  '.12) 

OtOO  2  FORM ATT  * <  <  OPEN  HE S S ACE > ) I  ERR •  '.14.'  AT  ',I2> 

0*01  3  FORHATT ' TTLOCF  HES SACE > > I ER R  -  '.14.'  AT  ' .  I  2  > 

0602  4  FORHATT ' <<CLOSE  HES SAGE >> I  ERR ■  '.14.'  AT  '.I2> 

0603  9  FORHATT ' <<WRITF  HE SSAGE > > I ER R ■  ',14.'  AT  ' ,  I  2  > 

0604  6  FORHATT ' <<READF  HES SAGE >> I  ERR ■  '.14.'  AT  ' , 1 2  > 

0609  10  FORHATT ' •••••TIHE  OUT  ERROR*****' 

0606  ♦/'••♦••SUFFER  OVERFLOW** **'.  I  6 > 

0607  END 

0  60  8  C* •**••**•••*•••*«•*«•**•*••*•****••*••*** •****•**•• . . 

0609  SUBROUTINE  TAPET  TAPELU ) ,  TO  CONTROL  A  TAPE  LU  BY  PROGRAM 

0610  C  :  A!t  A  GEN  PURPOSE  ROUTINE  TO  DO  HEAT  THINGS  TO  ANY  HAG-  C 
0611  C  NET  1C  TAPE  DRIVE.  C 

0612  C . - . --C 

0613  COHHON  LUTS) 

0614  INTEGER  T APELU . I BUF T 32 > 

0619  90  URITET LU,  100)  TAPELU 

0616  100  FORHATT' . ENTER  AN  OCTAL  FUNCTION  FOR  LU  ',12,'  1043' 

0617  ♦  /'  0000  •  CLEAR  THE  DEVICE- 

0618  ♦/'  0100  -  WRITE  AN  EOF' 

0619  ♦/'  0200  -  BACKSPACE  1  RECORD' 

0620  ♦/'  0300  >  FORWARD  SPACE  1  RECORD* 

0621  ♦  /'  0400  *  REWIHB' 

0622  ♦/'  0300  *  WRITE  AN  INTER-RECORD  GAP' 

0623  ♦/’  1300  -  FORWARD  1  FILE- 

0624  ♦/'  1400  >  BACK  1  FILE' 

0629  ♦/'  2700  ■  LOCATE  FILETCTU  ONLY)' 

0626  ♦  /'  7777  *  READ  i.  LIST  A  RECORD  OF  UP  TO  32  WORDS' 

0627  ♦/'  1111  «  TO  GET  OUT  OF  THIS  ROUTINE'/' - .') 

0628  READTLU. 101  )  IFC 

0629  101  FORMAT T  04  ) 

0630  IFT IFC.NE.1111B)  GO  TO  112 

0631  WRITETLU, 113) 

0632  113  FORHATT'  RETURNING  TO  HAIN  PROGRAM') 

0633  RETURN 

0634  112  IFT IFC. EO. 27008)  GO  TO  103 

0633  IFT IFC  EO .77778)  CO  TO  110 

0636  IC0DE>3 

0637  IOPl-0 

0633  103  CALL  EXECI I  CODE , TAPELU* I FC , l OP l  ) 

0639  GO  TO  30 

0640  103  WRITET LU. 106) 

0641  106  FORHATT 10X, 'ENTER  FILE  8  TO  LOCATE  1121') 

0642  READTLU. •  )I0P1 

0643  GO  TO  103 

0644  1  10  CALL  E HE CT1, TAPELU,  IBUF, 32) 

0  643  WRITET  LU.  1 1  1  >T  IBUFT  I  ), I -  1 , 32  ) 

0646  111  FORMAT  <  32 A2  1 

0647  GO  TO  50 

0648  EHD 


24 


APPENDIX  C  -  EXAMPLE: 

SIX  cm  METHANOL  PAN  FIRE 

REST! 

T  *0  0  003  IS  ON  CR  000  4 1  USING 

00036  BLKS  R-0000 

0  OO  1 

6  CH  PhN  FIRST  DATS  POIHT 

9002 

VERTICAL  IN) 

HOR  12  ON  T  Al  <«) 

SAHPLE  INT.  TINE 

0OO3 

14326880 

-  2 1  6 23 394 E- 02 

811007926-01 

0004 

31290054E-02 

43333414 E- 02 

. 46202880 

O00S 

7301 7776E-02 

57  00 1 27  2E- 02 

38174719 

0006 

-  22369782E-02 

. 40  437  38  0  E-03 

23592938 

0007 

- . 22369782E-02 

-  7 1 4  3373  7  E- 02 

31948799 

0009 

. 43054266E-02 

-  21625 39 4E-02 

69795833 

0  00  4 

-  91 599 92 9E- 03 

40437380E-03 

42270720 

0  010 

6372705  IE-01 

. 362491O6E-02 

. 41613333 

0011 

16927639 

. 40437 38 OE- 03 

. 82329601E-01 

0012 

16367457 

. 17079869E-02 

. 10239999 

0013 

1 7  66  7  62  9E -  02 

40457 380E- 03 

32767999 

0014 

3 1 29o  05  4E -  02 

. 435S5414E-02 

. 46202880 

0013 

3 1 290  05  4E -  02 

40437 38 OE- 03 

. 78233600E-0  1 

0016 

. 45054  266E -  02 

-. 21625394E-02 

.71434236 

0017 

. 41 84 962 OE -03 

. 404  37  380  E- 03 

.57671678 

0019 

. 3 1 290  05  4E -  02 

17079869E-02 

14581738 

0019 

. 3 1 290  05  4E -  02 

. 17079 869E- 02 

. 11796479 

0020 

3 1  290  05  4E -  02 

17079869E-02 

60293114 

0021 

. 3129005 4E -02 

40457 38 OE- 03 

.31293440 

0022 

. 26740372 

. 40437380E-03 

. 38830078 

0023 

. 4 1 84  9  620E- 03 

-  2 1 6 25394 E- 02 

. 64225280 

0024 

. 13424206 

-. 21625394E-02 

24739838 

0023 

37  633  680E- 0 1 

-. 88332607E-03 

. 3 1743996E-0  1 

0026 

43735323 

. 30249 106E-02 

44236797 

0027 

. 38962824E-02 

. 4335341 4E-02 

49479675 

0028 

7301 7776E-02 

. 302  49 10  6E-02 

.33084159 

0029 

22  369  782E -  02 

. 40437380E-03 

. 17121279 

0030 

44723734 

30249106E-02 

50133040 

0031 

87221 134E-02 

. 37001272E-02 

.41123837 

0032 

10  004017 

. 17079869E-02 

. 63487992E-0  1 

0033 

9  1733093 

. 302  49 10  6E-02 

. 10138080 

0034 

4  6847963 

. 40437380E-03 

. 18022400 

0033 

14  246948 

. 40437380E-03 

. 864235886-01 

0036 

13.434493 

- . 88332607E-03 

. 10813439 

0037 

. 31 290034E-02 

. 43555  4 1 4  E-02 

62386880 

0038 

6  3395241 

- . 88552607E-03 

43212794E-01 

0039 

4.7763882 

43535414E-02 

32767999 

0040 

6  6842194 

17079869E-02 

32767999 

0041 

7  6619654 

17079869E-02 

34886393E-01 

0042 

10.337299 

. 30249 106E-02 

. 66764802E-0  1 

0043 

. 13864873 

-  .  88 5 32 60 7 E- 03 

.  49807358 

0044 

. 4 1 849  620E- 0  3 

. 30249106E-02 

.  19496939 

0043 

31 290054E-02 

. 40437380E-03 

.  83558393 

0046 

. 15522179 

17079869E-02 

. 30135040 

0047 

. 55655  66  6E -  0 1 

-.591698486-02 

. 28671998 

0048 

. 7301777 6E -02 

. 30249106E-02 

. 50135040 

0049 

. 59640586E-01 

. 30249106E-02 

.  70431200 

0030 

41849620E-03 

- . 216255946-02 

. 42270720 

0051 

6167089 2E- 01 

170798696-02 

49807358 

0052 

4134962OE-03 

30249 106E-02 

. 69468153 

0033 

- . 91 599929E-03 

-.885326076-03 

.  65945596E-0  1 

0034 

. 58  962  824E- 02 

.  30 2 49 1 0 6 E-02 

.  13926399 

0033 

99823773 

. 404573806-03 

.31293440 

0056 

. 15692824 

-  88532607E-03 

.  55705595 

0057 

. 41849620E-03 

. 40437380E-03 

66846716 

0038 

. 43221126 

. 404373806-03 

. 80936936 

25 


0499 

.  39608161 

4  444 

. 17 66 762 9E -02 

4461 

12323640 

4462 

31 29403 4E-02 

4443 

3 1 29005  4£ -  02 

4464 

- . 2236978  2E -  02 

4463 

. 31 29 405 4E -02 

4446 

- . 73894039E-02 

4467 

- . 91 39992 9E -03 

4464 

. 44723734 

4469 

. 4 l B49620E- 03 

4474 

38 96282 4E -02 

4471 

17667629E-02 

4472 

17667629E-02 

4473 

. 37129  26  IE-01 

4474 

43034266E-02 

4473 

17667629E-02 

4476 

-  9 1 399 92 9E - 03 

4477 

. 4543426 6E -02 

4478 

. 63727031E-01 

4479 

41 84 962 0E -03 

4484 

. 3129005 4E -02 

4481 

17.288971 

4482 

3.4437102 

4483 

1.3831941 

4484 

45434266E-02 

4483 

. 72924316 

4436 

-  .  9 1 599  92  9E -  03 

4487 

16746795 

0488 

. 6167089 2E -01 

4489 

. 43454  26  6E -  02 

4494 

- . 22369782E-02 

4491 

. 4 1 849  62  OE- 03 

4492 

.  1 7  667  629E- 02 

4493 

. 1 7667629E- 02 

4494 

. 38962824E-02 

4493 

. 16927639 

4496 

19425631 

4097 

41849620E-03 

4098 

. 17667629E-02 

4499 

- . 91599929E-03 

4144 

. 12233069 

4141 

16367457 

4142 

.  31 Z90054E-02 

4143 

.  31294054E-02 

4144 

. 1 766  762  9E -  02 

4143 

4 1 84  9  620E -  03 

4146 

41849620E-03 

4147 

45<»54266E->2 

4143 

41849620E-03 

4109 

.  1 7  66  7  62  9E -  02 

4114 

.37641710 

4111 

.  17110023 

4112 

- . 22369782E-02 

4113 

. 6372705  IE -0 1 

4114 

. 35749697 

4113 

.  13120064 

4116 

-. 22369782E-02 

4117 

16567457 

4118 

.  4 1 949 62 OE- 03 

-  .  88332607E-03 

12779319 

. 4335541 4E-02 

54722559 

.  40437380E-03 

33587199 

30249  146E-02 

. 44468478 

.  37001272E-02 

38830078 

34249106E-42 

. 27033597E-01 

. 40457380E-03 

27361280 

17479 86 9E-02 

11960319 

17079 869E- 02 

74035672 

17079869E-02 

71434236 

. 30249146E-02 

. 33095676 

4335541 4E-02 

24412158 

37001272E-02 

61276162 

.  17  0  7986  9  E- 02 

. 54722559 

. 17079869E-02 

.  1  1468799 

57001272 E-02 

64552951 

. 17079869 E-02 

29982716 

40457380E-03 

63242233 

. 30249106E-02 

56688643 

. 435S5414E-02 

. 17940480 

. 17079869E-02 

. 72089601 

37001272E-02 

32604158 

-  883 32 60 7 E - 03 

14335999 

. 34249106E-02 

14663678 

. 34249 1 0  6E-02 

. 78643203E-0  1 

-. 883S2647E-03 

57343996 

. 40457380E-03 

. 50790393 

-  21625594 E-02 

31293440 

. 4335541 4E-02 

. 36864041 

. 302  49 106  E- 02 

.81592309 

17079869E-02 

57343996 

. 302  49 106E-02 

72089601 

.  40457380E-03 

. 86835191E-01 

. 17079869E-02 

. 1 1468799 

.  40437380E-03 

.21135360 

.  S7001272E-02 

.71106553 

. 70388645E-02 

. 37847036 

.  17079869E-02 

82903039 

3O2491O6E-02 

32358401E-01 

.  30  2 49 106  E- 02 

.31129599 

-  88552607E-03 

39157760 

.  40457380E-03 

. 66191351 

.  17079869E-02 

52428794 

.  3O2491O6E-02 

. 66355 1 99E-0 1 

.  17079869E-02 

.  73072636 

.  17079 86 9E-02 

.  23999680 

17O79869E-02 

4  i  779  1  96 

.  4M57380E-03 

.23756799 

17079969E-02 

. 38993919 

30249106E-02 

.41287678 

.  1 70  7986  9  E- 02 

83558393 

88552 60 7E-03 

.42926079 

-  .  88352607E-03 

68812799 

.  17079869E-02 

35717118 

. 30249 106E-02 

. 72089598E-0 1 

. 40457380E-03 

.  12533760 

40437380E-03 

22282240 

. 40457380E-03 

. 24248320 

.  17079 869E- 02 

. 16629758 

.  17079869E-02 

69468153 

0119 

. 36468273 

. 30249 106E-02 

. 50462723 

0  120 

82406759 

- . 88532607E-03 

12779519 

0  12  1 

3 1 290  05  4E -  02 

. 17079869E-02 

13892479 

0  122 

1621324? 

-  21625394E-02 

. 22609919 

0123 

-  22 369 782E - 02 

. 40437380E-03 

. 43008000E-0  1 

0124 

. 5964058  6E -  0 1 

1 70  79  86  9E- 02 

. 82247674 

0123 

17479461 

. 43555414E-02 

.13107198 

0126 

19629845 

57001 272E-02 

. 70778871 

012? 

31290054E-02 

17079 86 9E- 02 

. 40468478 

0  12  6 

6415981 1 

43335414E-02 

25339038 

0129 

. 31 29003 4£ -02 

70388 64 3E- 02 

. 73400319 

0130 

58962824E-02 

. 70588645E-02 

. 20479998 

0131 

40125716 

. 40437380E-03 

. 56360960 

'>132 

17 66 762 9E -02 

.  30 2 49 106 E- 02 

. 98303989E-0  1 

0133 

. 53 96282 4E -02 

-  88S52607E-03 

.  19824639 

•>134 

45054266E-02 

. 70588643E-02 

. 46202880 

'>135 

3 1 29»  05  4E -  02 

40  4  57  380  E- 03 

44564480 

>136 

59640386E -  01 

17O79869E-02 

27361280 

•>  l?7 

-  1334851 6E -01 

-  10731111 E- 01 

. 73728001 

>13  3 

4 1 34  9  62  oE- o3 

30249106E-02 

. 73728001 

>135 

3 1 290  054E -  02 

. 57O01272E-02 

. 17694718 

0140 

59640586E-01 

17079869E-02 

.64225280 

» 14  1 

-  22369782E-02 

30  249 106E- 02 

. 83558393 

0142 

62205803 

30  2  49 106  E- 02 

. 13400938 

0143 

40538666 

98532607E-03 

. 78643191 

0144 

3 1 290  054E -  02 

.  30 2  49  10 6 E- 02 

. 44364480 

0145 

45054  266E-02 

1 70  7986  9  E- 02 

. 70451200 

0146 

4184962 OE -03 

.  70588645E-02 

. 43581438 

0  14? 

4184962 OE- 03 

43335  41 4  E- 02 

80936956 

0148 

4  t  84  9  620E -  03 

. 40457380E-03 

. 47513598 

0149 

. 3129OO54E-02 

. 30249106E-02 

. 93388796E-0 1 

0150 

. 87221 134E-02 

-  885526O7E-03 

. 59310079 

0131 

3 1 29O054E- 02 

4353341 4E-02 

80936956 

0132 

58962824E-0 2 

30249 106E-02 

. 60293114 

0153 

-  9 1 39  9  92  9E- 03 

40437380E-03 

. 36208636 

0134 

1 766  7  629E- 02 

. 40457380E-03 

.  13646720 

0133 

31290054E-02 

17079869E-02 

. 1 1304960 

0136 

41849620E-03 

- . 883326O7E-03 

39976960 

0137 

17667629E-02 

. 57001272E-02 

36864001 

0138 

44725734 

. 37001 272E-02 

. 54067194 

0159 

7301 7776E-02 

. 57001272E-02 

.71434236 

0  160 

.  31 290034E-02 

. 4333541 4E-02 

. 27033397E-01 

0161 

12235069 

17O79869E-02 

13352960 

0162 

3 1 290054E -  02 

. 30249 106E- 02 

. 77004790E-0  1 

0163 

65809533E-01 

30249106E-02 

50135040 

'>164 

15322179 

.  1707986  9E-02 

.52101111 

0165 

-  9 1 399  92  9E -  03 

885326O7E-03 

. 78315520 

0166 

.  4 l 849 62 OE - 03 

- . 88352607E-03 

. 28508157 

016? 

73017776E-02 

. 43535414E-02 

. 73400319 

0163 

4305  4  26  6E -  02 

.  40 437 38 OE- 03 

. 77332473 

0  169 

-  9 1  39  9  92  9E  -  o3 

- . 885526O7E-03 

. 10158080 

0170 

-  9 1 599  92  9E -  03 

.  40437380E-03 

. 82903039 

0171 

. 84763923 

40437380E-03 

72499201 E-0  1 

0172 

38448751 

.  70 3 8864 5E- 02 

47513598 

0173 

4194962 OE -03 

.  4O437380E-03 

.  93027 1 93E-0  1 

0174 

. 59640386E-01 

.  40437380E-03 

.  45547513 

0173 

7301 7776E-02 

.  43535414E-02 

.  49151999 

0176 

.  61  670892E-01 

170  79  86  9E- 02 

. 33382401E-01 

0177 

41 84962 0E-03 

40437380E-03 

26369756 

0178 

6620221 1 

30249 106E-02 

56360960 

0179 

4184962 o£-«3 

88552607E-03 

31293440 

0184 

61670 89 2E -Oi 

40437380E-03 

. 43875198 

0181 

91138331 

40  4  57  38  o  E- 0  3 

47513598 

0192 

3 1 290  05  4E -  02 

40457380E-03 

33259517 

0183 

17110023 

43553414E-02 

24412158 

0184 

. 41 849620E-03 

43555414E-02 

21954557 

0183 

45034  26  6E -  02 

37001272E-02 

. 66191351 

0186 

17667629E-02 

40  4  37  38  0  E-  03 

16957438 

0187 

-  9 1 399 329E- 03 

40  437380  E*  03 

53247996E-0 1 

0189 

4184962OE-03 

17079869E-02 

212991 98 E-01 

0  139 

17667629E-02 

88  5  52  60  7  E-  03 

. 67502081 

0  190 

45054  26  6E -02 

40437380E-03 

. 80936956 

0191 

41 849620E-03 

43333414 E- 02 

22609919 

0192 

-  61 20 602 6E -02 

37001272E-02 

34078717 

0193 

- . 2236  978  2E -  02 

30  2  49 106  E- 02 

14499840 

0194 

-  1354831 6£ -01 

30249106E-02 

69468153 

0193 

- . 91 399929E-03 

40437380E-03 

. 39321595 

0196 

28436261 

302  49 106  E- 02 

64307 1 98E-0  1 

0197 

43034  266E -  02 

4335341 4E-02 

. 33751041 

0198 

- . 9 1  599929E -  03 

40437380E-03 

. 11960319 

0199 

. 17667629E-02 

17079869E-02 

23592958 

0200 

39  64  0  38  6E -  0 1 

17079869E-02 

.69140470 

0201 

. 41849620E-03 

40457380E-03 

. 901 12001E-01 

0202 

. 27013636 

30249 106E-02 

. 78315520 

0203 

. 17667629E-02 

88332607E-03 

. 12779519 

0204 

. 9804461 9E-01 

2 1 6 23 394  E- 02 

. 79953921 

0203 

39640386E-01 

17079869E-02 

79953921 

0206 

0207 

0208 

0209 

0210 

0211 

0212 

0213 

0214 

0213 

0216 

0217 

31290034E-02  . 30249 106 E- 02 

UNWEIGHTED  MEAN  VELOCITIES* 

<  CNTR  1>  .33307307  K/S 

< CNTR  2 >  .  1 793 1 876E-02  N/S 

RESULTANT  VECTOR* 

MAGNITUDE  •  0.33307605 

INCLUDED  ANCLE  *  0.19294963  FROM  VERTICAL 

. 72908789E-0 1 
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